Modified BCS mechanism of Cooper pair formation in narrow energy bands of special 

symmetry 
I. Band structure of niobium 
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The superconductor niobium possesses a narrow, roughly half-filled energy band with Bloch func- 
tions which can be unitarily transformed into optimally localized spin dependent Wannier functions 
belonging to a double-valued representation of the space group 0\ of Nb. The special symmetry 
of this "superconducting band" can be interpreted within a nonadiabatic extension of the Heisen- 
berg model of magnetism. While the original Heisenberg model assumes that there is exactly one 
electron at each atom, the nonadiabatic model postulates that the Coulomb repulsion energy in nar- 
row, partly filled energy bands is minimum when the balance between the bandlike and atomiclike 
behavior is shifted as far as possible towards the atomiclike behavior. Within this nonadiabatic 
Heisenberg model, the electrons of the superconducting band form Cooper pairs at zero tempera- 
ture. Just as in the BCS theory of superconductivity, this formation of Cooper pairs is mediated 
by phonons. However, there is an important difference: within the nonadiabatic Heisenberg model, 
the electrons in a narrow superconducting band are constrained to form Cooper pairs because the 
conservation of spin angular momentum would be violated in any normal conducting state. There is 
great evidence that these constraining forces are responsible for superconducting eigenstates. That 
means that an attractive electron-electron interaction alone is not able to produce stable Cooper 
pairs. In addition, the constraining forces established within the nonadiabatic Heisenberg model 
must exist in a superconductor. 

PACS numbers: 74.10.+V, 74.20.-z, 74.20. Fg, 74.25. Jb 



I. INTRODUCTION 



The original Heisenberg model of magnetismEl is de- 
fined by the assumption that there is exactly one electron 
at each atom of a metal. Within the nonadiabatic Heisen- 
berg model (NHM) as proposed in previous papersHa the 
great success of the Heisenberg model is interpreted by 
introducing three new postulates given in Sec. 0. 

As a consequence of the postulates of the NHM, the 
ground state of the electron system consists of config- 
urations in which as many as possible atoms are occu- 
pied by exactly one electron. However, this state with 
the highest possible atomiclike character cannot be de- 
scribed within the adiabatic (or Born-Oppenheimer) ap- 
proximation. Therefore, the adiabatic localized states 
represented by Wannier functions (Wfs) are replaced by 
nonadiabatic localized states depending on an additional 
quantum number v. This new quantum number labels 
different states of that part of the motion of the center of 
mass of a localized state which nonadiabatically follows 
the motion of the electron. 

Hence, within the NHM the electronic motion is cou- 
pled to the nonadiabatic motion of the centers of mass, 
and, via this coupling, the electrons couple to each other 
and to the phonons. The character of the resulting 
electron-electron, electron-phonon, or spin-phonon inter- 
action depends crucially on the symmetry properties of 
the nonadiabatic localized functions. In this paper, a new 
type of spin-phonon interaction will be derived within the 



NHM. 

According to the third postulate of the NHM, the 
nonadiabatic localized functions have the same symme- 
try as the best localized Wfs related to the narrowest 
partly filled energy bands of the metal under consider- 
ation. Thus, any application of the NHM starts with 
a group-theoretical examination of the symmetry of the 
best localized Wfs which is determined by the symme- 
try of the Bloch inactions in the band structure of the 
considered metal. Em 

Paramagnetic iron possesses a narrow, roughly half- 
filled "ferromagnetic" band with Bloch functions which 
can be unitarily transformed into optimally localized Wfs 
belonging to a one-dimensional corepresentation of the 
magnetic group M — IA/mm'm' of the ferromagnetic 
state. In Ref. [|it was shown that if the postulates of the 
NHM are satisfied within the ferromagnetic band, then 
the electron spins form a spin structure with the mag- 
netic group M. n 

In an earlier papers I have shown that chromium pos- 
sesses a narrow, roughly half-filled "antiferromagnetic 
band" with Bloch functions which can be unitarily trans- 
formed into Wfs that are symmetry-adapted to the mag- 
netic group M = PiA/mnc of the spin-density wave state 
existing in this metal. Again it can be shown that «/the 
postulates of the NHM are satisfied within the antiferro- 
magnetic band, then the electron spins form a spin struc- 
ture with the magnetic group M of the antiferromagnetic 
state.B 

In this paper I consider the band structure of the su- 
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perconductor niobium as depicted in Fig. El in particular 
the band denoted by the dotted line. This band is char- 
acterized by the four representations 



T' 25 , H 25 , N 2 , and P 4 



(1.1) 



of the space group 0\. The Bloch functions of this 
single band cannot be unitarily transformed into usual 
Wfs which both are best localized and belong to a rep- 
resentation fir-Ph Si 
dimensional BQ 



since r' 25 , H 25 , and P4 are not one- 




FIG. 1. Band structure of Nb after Mattheis (Ref. . The 
dotted line denotes the superconducting band. 

However, such a unitary transformation becomes pos- 
sible in a natural way when we allow that the Bloch 
states of the considered band have k dependent spin 
directions. With this generalization we may construct 
from the Bloch functions of a band with the representa- 
tions (1.1) spin dependent Wfs which belong to a double- 
valued representation of Oh] see appendix 0. 

Consider a metal with the (bec) space group 0\ and 
one atom per unit cell. An energy band of this metal 
is called "superconducting band" (cr band) if the Bloch 
functions of this band can be unita rily t ransformed into 
spin dependent Wfs [as given in Eq. ( A14)] which ar e bes t 



locali zed, symmetry-adapted [according to Eqs. ( A15| ) 



and ( |A16| )1 to the paramagnetic space group M p and 
centered at the atoms. All the a bands in bec metals are 
given in Table fl. According to this definition, the consid- 
ered band ( fO ) of Nb is a a band of the symmetry type 
of band 2 in Table |; see appendix |A[ 

In Sec. H the three postulates defining the NHM will 
be given and in the following section [ill] I shall show that 
these postulates yield a special spin-phonon interaction 
in narrow, partly filled a bands. At zero temperature, 
this spin-phonon interaction constrains the electrons of 
the a band in a new way to form Cooper pairs because 
the conservation of spin angular momentum would be vi- 
olated in any normal conducting state. That means, if 
the postulates of the NHM are satisfied within a a band, 
then the electrons of this band form Cooper pairs at zero 
temperature. 

The mechanism of Cooper pair formation within the 
NHM is nearly identical to the familiar mechanism 



presented within the Bardeen-Cooper-Schrieffer (BCS) 
theory.E3 The important difference between both mecha- 
nisms will be given in Sec. IV. 

In appendix |a| the spin dependent Wfs related to a a 
band will be defined and the essential properties of the 
related nonadiabatic localized functions will be given. 



II. NONADIABATIC HEISENBERG MODEL 

In this section the equations defining the NHM will 
be given together with a short substantiation of the new 
model. For a more detailed describtion see Sec. II of 
Ref. H 

Consider a single partly filled energy band in a metal 
with one atom in the unit cell. In this section, this band 
need not yet be a cr band. Let 

H = H H F+H cb (2.1) 
be the electronic Hamiltonian of this band with H^p and 



H C b = , ^{fi,m 1 ]f 2 ,m 2 \Hcb\T[,m' 1 ]f! 1 ,m' 2 ) 



E 

f,m 



(2 2) 

"rimrT 2 m 2 " J 2 m 2"' J i m 'i V ' ' 

representing the Hartree-Fock and Coulomb energy, re- 
spectively. The fermion operators ct, and Cf m create 

and annihilate electrons in the localized states \T, m) rep- 
resented by (spin dependent) Wfs w m (r — f, t) (with the 
spin coordinate t) which are assumed to be situated at 
the atoms and to form a complete basis within the con- 
sidered band. Other contributions to H from the elec- 
trons not belonging to this band are neglected even as 
are spin-orbit effects. 
Hcb may be written as 

Hcb = H c + H ex + H z , (2.3) 

with the operator of Coulomb repulsion H c containing all 
the matrix elements of Hcb with 

fx = f[ and f 2 = % (2.4) 

the exchange operator H ex containing the matrix ele- 
ments with 

fx = f 2 and f 2 = f{, (2.5) 
and H z comprising the remaining matrix elements. 
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TABLE I. Double-valued representations ii- of the four 
superconducting bands in the bec structure with notations of 
Ref. El. 
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Unlike H c or H ex , the operator H z generates virtual 
transitions between adjacent localized states which in 
narrow bands shift the balance between the bandlike and 
atomiclike behavior toward the bandlikc character. The 
NHM assumes that these transitions are energetically 
unfavorable in sufficiently narrow, partly filled energy 
bands. Hence, the first postulate of the NHM states that 
the relation 



(G\H\G) > (G'\H'\G') 



(2.6) 



is satisfied in the narrowest half-filled energy bands of 
the metals. The new Hamiltonian 



H' = Hhf + H c + H e 



(2.7) 



is obtained from H by putting H z — 0, and \G) and |G") 
stand for the exact ground states of H and H', respec- 
tively. 

The particular form of the matrix elements of H z shows 
that it represents a short-ranged interaction which cru- 
cially depends on the exact form of the localized func- 
tions. This fact suggests that only small changes of 
the localized electronic orbitals are required to prevent 
the transitions generated by H z . However, such mod- 
ifications of the electronic orbitals do not exist within 
the adiabatic approximation because these modifications 
yield charge distributions within the localized states be- 
ing symmetric with respect to the lattice on the average 
of time, but not at any moment. Consequently, the pos- 
itive ions become permanently accelerated in varying di- 
rections. Therefore, within the nonadiabatic model, we 
replace the (adiabatic) localized states represented by the 
Wannier functions by more realistic nonadiabatic local- 
ized states, 



\T,m,p), 



(2- 



which take into account the motion of the nuclei. The 
new quantum number v labels different states of motion 
of the centre of mass of the nucleus and the electron oc- 
cupying the state \T, m, v). 

The nonadiabatic Hamiltonian H n may be written as 



H n = H 



HF 



H Cb 



(2.9) 



where the Coulomb interaction now has the form 



H cb = y"^(7i) mi,n; T 2 , m 2 , n\Hcb\T[, m'^n; 7%, m' 2) n) 

(2.10) 



nf nf n n 
f imi °f 2 m 2 °fi m' 2 C f \ m[ ' 



The new fermion operators c~ 



and cS create and an- 
nihilatc electrons with crystal spin m in the nonadiabatic 

n labels the nonadia- 



localized states |T, m, n) and v 
batic states which satisfy the following Eq ( p. 12 ). The 
new coordinate q represents that part of the motion of 
the nuclei which nonadiabatically follows the motion of 
the electron occupying the state \T, m, n). 



The second postulate of the NHM states that the tran- 
sitions generated by H z are artifacts of the adiabatic 
approximation and do not ha ppen in the (true) nona- 
diabatic system if relation (2.6) is satisfied. That means, 
the NHM assumes that there exist nonadiabatic localized 
functions 

(f,t,q\T,m,n) (2.11) 

satisfying 

(T 1 ,mi,n;T 2 ,m2,n\Hcb\T{,m[,n;T2,m' 2 ,n) =0 (2.12) 
if 

{fi.fa}^^,^}. 

According to the third (and last) postulate, these nona- 
diabatic localized functions have the same symmetry as 
the Wfs w m (f-f,t). 



III. APPLICATION OF THE NONADIABATIC 
HEISENBERG MODEL TO 
SUPERCONDUCTING BANDS 

A. The Coulomb interaction in a narrow a band 

Now assume the narrow, roughly half-filled energy 
band considered in the preceding section to be a a band. 
That means that now the Wannier functions w m (f—T, t) 
in Eq. (O) are given by Eq. (|A14[) . 



We first examine the commutation properties of the 
operator H 1 [given in Eq (2.7)]. [It should be noted that 
the symmetry of H' depends on the symmetry of the 
Wannier functions w m {f—T, t), whereas the symmetry of 
the complete Hamiltonian H is independent of the sym- 
metry crfthe use d ba sis functions.] As a consequence of 
Eqs. ( A25 ) and ( A26 ) it follows that H 1 commutes with 
the symmetry operators P{a) of the space group G, 



[H',P(a)] = for a e G, 



(3.1) 



and with the operator K of time inversion, 

[H',K]=0. (3.2) 



In Eq. (All) the symmetry operators S(a) of the elec- 
tron spin are introduced. As shown in appendix the 
Coulomb interaction 



H'cb = Hc + He. 



(3.3) 



belonging to H' does not commute with the operator 
S(a), 



[H' cbl S(a)}^0 
for at least one a S Gq-Q 



(3.4) 
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Now assume the three postulates of the NHM to be 
satisfied in the a band of niobium. According to the 
third postulate, the nonadiabatic Hamiltonian H n has 
the same commutation properties as H' . Hence, H n com- 
mutes with the operators P{a) and K, 



and 



[H n ,P(a)] =0 for aEG 



[H n , K] = 0, 



(3.5) 



(3.6) 



where the action of the oper ators P(a) in the nonadia- 
batic system is given in Eq. ( A41 ). 

As a consequence of Eq. ( |3.4| ) also t he no nadiabatic 
Coulomb interaction ifg 6 [given in Eq. ( |2.1C| )1 does not 
commute with all the operators S(a), 

[HZ b ,S(a)}^0 (3.7) 

for at least one a E G . 

The first two equations (^^) and (3.6) show that H r ' 



has the correct symmetry of t he p aramagnetic group M p 
of niobium. Especially, Eq. ( |3.6|) indicates that H n has 
a nonmagnetic groun d sta te. 

The last equation (3.7) demonstrates that H n repre- 
sents a system in which the electron spin angular mo- 
mentum is not a conserved quantity. Within the nona- 
diabatic system, we need not require that the electron 
spin is conserved because the nonadiabatic fermion op- 
erators c™ and d% are no longer labeled by the spin 

Tm Tm ° . . 

quantum number s. However, even in the nonadiabatic 
system there should exist a conserved quantity related to 
the conservation law of angular momentum. Thus, the 
equations 



[H,S(a)} = for a E 0(3), 



(3.8) 



which hold in the adiabatic system, should be replaced 
by analogous equations 



[H n ,M{a)\ =0 for a e Go 



(3.9) 



in the nonadiabatic system. 

The operators M(a) are defined in Ref. |[ They act on 
the quantum number m of the nonadiabatic fermion op- 
erators in the same manner as the operators S(a) act 
on the spin quantum number s of the adiabatic fermion 
operators ct. . Therefore, the operators M(a) may be 
called the "operators of crystal spin" and m may be called 



the quantum number of the crystal spin. This is in anol- 
ogy to the wave vector k of the Bloch functions which is 
sometimes referred to as "crystal momentum" in order 
to distinguish it from the momentum p. 

The operator H' cb has matrix elements with 



mi + 777,2 7^ Till + m 2 



(3.10) 



because the matrix [f sm {k)\ in Eq. ( |A9| ) is not indepen- 
dent of k in a a band. Therefore^, H' Cb does not conserve 
the crystal spin. According to the third postulate of the 
NHM, also the nonadiabatic Coulomb interaction W c 
a narrow a band does not conserve the crystal spin, 



l Cb 



[HZ b ,M(a)}^0 



(3.11) 



for at least one a E G . 

In order to interprete this result remember that in the 
(true) nonadiabatic system the electrons no longer move 
on rigid orbitals in the average potential of the other elec- 
trons, but move in a potential depending on which of the 
adjacent localized states are occupied and on the present 
positions of these electrons. These modified electronic 
orbitals yield charge distributions within the localized 
states being symmetric with respect to the lattice on the 
average of time, but not at any moment. As already 
stated above, the positive ions become permanently ac- 
celerated in varying directions within this fluctuating po- 
tential of the electrons. 



Equation (3.11) indicates that in a narrow a band 
the positive-ion lattice is accelerated in such a way that 
phonons are excited (or absorbed) . This interpretation is 
corroborated by the fact that the acoustic phonons in a 
cubic crystal may be unitarily transformed into localized 
phonon states |T, 1} which transform according to 

M(a)\T,l) =^2Din(a)\f,l') for a € G (3.12) 

by application of the operators M(a) of the crystal spin.i 
The matrices [Din(a)] belong to the three^iimcnsional 
representation Tis in every cubic crystal.E3 Therefore, 
the localized phonons |T, I) carry the crystal-spin mo- 
mentum 1 • h and the quantum number I may be identified 
as crystal-spin quantum number I = —1, 0, +1. 

T here is further evidence for this interpretation of 
Eq (3.11) because by the mere addition of symmetrized 
boson operators we may construct from H^ b an interac- 
tion 



HZl = ^(f 1 ,; 1 ;f 2 ,; 2 ;f 1 ,m 1 ,n;f 2 ,m 2 ,n|iJc b |T 1 ,mi,n;^,m' 2 ,n)^ bt c" + & c n f c n f + H.c. (3.13) 
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which conserves the crystal spin, 

[H%l,M(a)] = for a <= G . 



(3.14) 



The boson operators bL, and brp, create and annihilate 
it 



localized phonons \T,l) with crystal spin / = — 1,0, +1 



at the position Tl^The complicated prove of Eq. ( 3.14 ) 
is given in Ref. 0.11-3 

Though conserves the crystal-spin angular mo- 

mentum, it still does not conserve the electron spin (since 
the nonadiabatic fermion operators are not labeled by the 
spin quantum number s). Hence, we still have 



[H%,S(a)]?0 



(3.15) 



for at least one a £ Go- This equation indicates that it is 
the electron spin which couples to the phonons. Hence, 
the electron-phonon coupling produced by H^l may be 
called "spin-phonon interaction" . 

The interaction HqI replaces the Coulomb interaction 
Hqi, in a narow a band. The complete nonadiabatic 
Hamiltonian now may be written as 

H n = H HF +H^ b +H ph , (3.16) 

with H p h denoting the operator of phonon energy. 



B. Formation of Cooper pairs 



The result of the preceding section [III A| may be sum- 
marized by an if-then statement: if the postulates of the 
NHM are satisfied in a narrow, partly filled a band, then 
the electrons are coupled to the phonons by the nonadi- 
abatic Coulomb interaction Hq^. In this section I shall 
show that the new interaction necessarily produces 
Cooper pairs at zero temperature because in a normal 
conducting state the conservation of angular momentum 
would be violated. 

Since Hgifc depends on phonon operators, a certain 
number of crystal-spin- 1 phonons is excited in the ground 
state of the nonadiabatic Hamiltonian H n at any temper- 
ature. At zero temperature we may assume that these 
phonons are virtually excited, i.e., each phonon pair is 
reabsorbed immediatedly after its generation. Thus, at 
zero temperature it should be possible to approximate 
the operator H n given in Eq. (3.16) by a purely elec- 
tronic operator 



H° = H 



HF 



H Cb 



(3.17) 



using a canonical transformation-ip analogy to the proce- 
dure first discussed by Fr6hlich.Ej Now the Hamiltonian 
no longer depends on the boson operators bt, and bf v 

Just as HqI, the operator H^, b conserves the crystal 
spin, 



[H% h ,M{a)) = for a e Go- 



Thus, the ground state |G°) of H° satisfies (in the sim- 
plest case) the equation 



M(a)\G°) = \G°) for a e G . 



(3.19) 



H° does not depend on boson operators. Conse- 
quently, phonons are not excited within \G ) though 
also at zero temperature the electron spins couple to 
the motion of the centers of mass and, hence, to the 
phonons. However, the phonon system does not store 
angular momentum but mediates (via virtual phonons) 
a new electron- electron interaction which now is repre- 
sented by Hg h . Thus, the electron system becomes adia- 
batic when H n is replaced by H°. The localized electron 
states now are represented by the adiabatic spin depen- 



dent Wfs (r,t\T,m) given in Eq. ( |A14| ) and the electron 
spin is a conserved quantity, 



[H° b ,S(a)}=0 foraeO(3), 

yielding 

S(a)\G°) = \G°) for a € 0(3). 
Equation ( A20| ) may be written as 



4 

km 



BvK 

E 



(3.20) 



(3.21) 



(3.22) 



and Eq. (A9) yields 



(3.23) 



where the (adiabatic) fermion operators ct and ct cre- 

km ks 

ate Bloch electrons with crystal spin m and spin s, re- 
spectively. 

We now show that both equations ( 3. 19] ) and (3.21) are 
satisfied if and only if the electrons form Cooper pairs. 
Consider states of the form 



where the new operators 

dt = ct c 1 - - 

r k fe-ff -ki}. 



•4 |o>, 



ct c t . 

fe^ -kit 



(3.24) 



(3.25) 



create symmetrized Cooper pairs. The crystal-spin quan- 
tum numbers m — and m — — ^ are denoted by dou- 
ble arrows ff and Jj-j respectively. L denotes the number 
of electrons. 

The operators 0t are basis functions of the identity 
representation I\ of Gq. Thus, we have 



M{a)0t = f3t for a e G , 



(3.26) 



(3.18) and hence 
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M{a)\Cp) = \Cp) foraeGo- 



(3.27) 



The operators /?! may be transformed into the s rep- 



resentation. With Eq. (3.23) we obtain 

4 = (av*)4 t + AV*)i 



(3.28) 



where the spin quantum numbers s — + ^ and s = — i 
are denoted by single arrows j and J,, respectively. 
This equation ( |3.28| ) leads with Eq. ( |A33| ) to 



#1 = c t c f - - ct c f - . 

A fe fe| -fcj. fcj. -fct 



(3.29) 



Thus, within the s representation the operators fit form 
basis functions of the identity representation of 0(3). 
From this important result it follows that, in addition 
to Eqs. Q3.26D and ( p^27| ), we have 



S(a)(3t = (3l for a G 0(3) 



(3.30) 



and 



which represent scattering processes changing the wave 
vectors k of the Bloch states and which commutes with 
both S(a) [for a G 0(3)] and M(a) [for a G G ] are 
given by 



since the equations (3.26) and (3.30) may be written as 



and 



[fit, M(a)} =0 for a G Go 



[/?t 5(a)] = foraeO(3). 



(3.33) 



(3.34) 



Th us, as a consequence of the equations ( 3.18| ) 
and ( 3.20 ), H^ b has the form 



H c b = J2&\H Cb \k')f3l(3p. 



(3.35) 



k,k' 



Hence, in the system represented by Hq the Coulomb 
interaction H^ b is strongly k and s dependent since 



(fci, si;fc 2 , S2\H% b \k[, s[;k 2 , s 2 ) = 



(3.36) 



S(a)\Cp) = \Cp) for a G 0(3). 



(3.31) 



Consequently, if the ground state |G°) of H° is a linear 
combination of the states \Cp), i.e., if all the elec tron s 
in | G°) fo rm Cooper pairs, then both equations ( 3.1£ ) 
and (3.21) are satisfied. 

The sec ond i mpor tant result is that the equa- 
tions (3.19) and (3.21) cannot be satisfied at the same 
time when the electrons in |G°) do not form Cooper 
pairs because the matrix [f sm (k)] in Eq. ( |A9| ) cannot be 
chosen ind epen dent of k in a a band; see appendix [A|. 
Equation (A33) is a consequence of the time- inversion 
symmetry of th e ele ctron-ph onon system. Hence, the 
step from Eq. (|3.28|) to Eq. (3.2S) is possible only for 



the symmetrized Cooper pairs (3l which consist always 
of both, an operator ct and its time-inverted operator 

±ct B • 

— k,—m 

An essen tial p ropert y of H% h can be derived from the 
equations ( 3.18| ) and fl3.20p . Assume that H^ b may be 
represented by a two-electron interaction, 



H Cb 



E 

k,s 

xct 



(fci, si; k 2 , S2\H^ b \k[, s[; k' 2 , s' 2 ) 



ct 



k 2 s 2 k 1 s 1 



(3.32) 

The only fermion operator combinations of the form 



E 



d(sx, s 2 ; s[, s 2 )ct 



ct 



fclSi k 2 S2 k' 2 



for fci 7^ —k2,k[ ^ —k' 2 ,si ^ —s 2 , or s'i ^ —s' 2 . (In ad- 
dition, Hg b has nonvanishing matrix elements belonging 
to scattering processes not changing the wave vector k. 
Here, these diagonal and exchange matrix elements are 
not considered.) 



C. Calculation of the transition temperature T c 

Within the NHM, the superconducting transition tem- 
perature T c may be calculated in the framework the BCS 
theory when, first, the nonadiabatic operator H^ b of 
Coulomb interaction is approximated by the adiabatic 
operator of Coulomb interaction, Hcb, and the usual 
spin-independent operator of electron-phonon interac- 
tion, H e _ p h, 



(3.37) 



with 



H e _ ph = J2(k, k\H e - ph \k[)bict i>s cl [ s + H.c, (3.38) 



and when, second, the approximations of the BCS theory 
are applied Hcb + Hc-ph + H p h in order that the familiar 
equation 



T r = 1.14-0 .e-W v 



(3.39) 







may be derived. N(Ep), V and 9 are the density of 
states of the electrons at the Fermi level, the effective 
electron-phonon interaction, and the Debye temperature, 
respectively. 

There is no definite indication of the approxima- 



tion (3.37) leading to a markedly false calculation of the 
free energy of the ground state and, hence, of T c . It is 
true that H^f h represents a two-electron two-phonon pro- 
cess. This, however, does not mean that the two phonons 
become excited (or absorbed) at the same time. The spe- 
cial form of H£1 should be interpreted by stating that, 
after a one-electron one-phonon process, the probability 
for another correlated one-electron one-phonon process 
rises. The accurate (two-electron, two-phono n) for m of 
HqI serves only to show the valid ity o f Eq. (3.14) and 
hence to deri ve th e important Eq. ( 3.36 ). 

Eq uatio n ( 3.36 ) has no influence on the BCS equa- 
tion (3.3E) because the matrix elements of H^, b which 
vanish in the NHM are argued away in the BCS theory, 
too. Hence the only contribution of the NHM to the 
BCS equation is that N(Ep) now stands for the density 
of states N°{Ef) of the electrons belonging to the jth a 
band. Within the NHM equation ( 3.39 ) reads as 



Tl = 1.14 ■ 9 ■ e -V"f 



(3.40) 



The NHM thus modifies the interpretation of the BCS 
formula in that now the parameter N(Ef) is only the 
partial density of states NJ{Ef) of the jth a band. 

The index j is introduced since the considered metal 
may possess more than one a band in its band structure. 
In this case, Tj! i s calc ulated for e ach a band seperately, 
since equations ( 3.19| ) and (3.21) are satisfied within a 
a band only. Hence, each of the a bands forms its own 
superconducting system with its own transition temper- 
ature. 



IV. DISCUSSION 

The result of this paper may be summarized by an if- 
then statement: if the three postulates of the NHM are 
satisfied within a narrow a band, then the electrons form 
Cooper pairs at zero temperatur. Just as in the BCS the- 
ory, the coupling of the electrons is mediated by phonons. 
However, the formation of Cooper pairs is constrained in 
a new way by the conservation law of crystal-spin angu- 
lar momentum: within a narrow, partly filled a band the 
crystal-spin angular momentum is conserved at zero tem- 
perature if and only if the electrons form Cooper pairs. 

From a mathematical point of view, the difference be- 
tween the BCS mechanism of Cooper pair formatio n an d 
the mechanism within the NHM is seen from Eq ( 3.36 ). 
In both the BCS t heory and the NHM the matrix ele- 
ments given in Eq ( |3.36 ) may be disregarded. However, 
the methods to argue away these matrix elements are 
quite different. 



In the NHM, the matrix elements in Eq. (3.36) may 
be clearly disregarded because they vanish. Within the 
BCS theory, on the other hand, the corresponding matrix 
elements of the the effective electron-electron interaction 

V = y^(fci, si, k 2 , s 2 \V\k[,s 1 ; k' 2 ,s 2 ) 

k . s 

xct ct cr, cr, (4.1) 

fclSl k 2 S 2 k 2 S2 k l Sl 

(dcrivedO from the spin-independent electron-phonon in- 
teraction) do not vanish from the beginning. But after 
the introduction of the familiar BCS wave function \<p) 
they may be disregarded because they vanish in the space 
spanned by the BCS function, 

(h,s 1 ;k 2 ,s 2 \c!))(<t>\V\ ( l))(<l>\k[,s[;k 2 ,s 2 ) = (4.2) 

for fci 7^ — k 2 ,k[ 7^ — k 2 ,si 7^ — s 2 , or s[ ^ ~s 2 . Conse- 
quently, these matrix elements may be disregarded if and 
only if the BCS Hamiltonian 



H 



BCS 



H 



HF 



V 



(4.3) 



possesses eigenstates in that part of the Hilbert space 
which is spanned by the BCS wave function. Hence, 
the BCS theory uses the tacit assumption that Hbcs 
possesses eigenstates in which the electrons form Cooper 
pairs. 



The operator H° [given in Eq (3.17)] clearly has eigen- 
states in which the electrons form Cooper pairs because 
its complete interaction term H^ b satisfies Eq. (3.36). 
However, it is conceivable that Hbcs does not possess 
such eigenstates since H° and Hbcs act in different sub- 
spaces V° and Vbcs of the Hilbert space. 

The BCS operator V acts in a space Vbcs in which 
each Bloch electron is described by the four quantum 
numbers k x ,k y , k z and s. In the space V°, on the other 
hand, each electron p air i s described by four quantum 
numbers because Eq. ( 3.36| ) is valid. Thus, V° and Vbcs 
essentially differ because the number of degrees of free- 
dom of the electronic motion in V° is one-half that in 
Vbcs- 

In a classical system, the number of degrees of free- 
dom of any state of motion of N particles is given by the 
number of independent coordinates on which the Hamil- 
tonian function depends, i.e., it is prescribed. If this is 
also true in a quantum system, then Hbcs has only un- 
paire d eig enstates and it is the operator H as given in 
Eq. (3.17) that has eigenstates in w hich the electrons 
form Cooper pairs. Equation ( 3.36|) is an equation of 
constraint which may be interpreted in terms of "con- 
straining forces!!, which constrain the electrons to form 
Cooper pairs Bll3 

Hence, it cannot be excluded that any electron-boson 
interaction which produces stable Cooper pairs within 
the BCS limit has two characteristic features: 

(1) It is attractive in the sence that is yields an energy 
minimum in the weak coupling (BCS) limit and 
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(2) the complete interaction term H^ b (involving both 
the Coulomb intera ction and the electron-boson interac- 
tion) satisfies Eq. (3.36) in order that the Hamiltonian 



with 



possesses superconducting eigenstates. 

In other words, there is evidence that the BCS theory 
yields the absolute energy minimum in the Hilbert space 
only if the formation of Cooper pairs is additionally con- 
strained by the constraining forces existing in a narrow, 
partly filled a band. 

As a consequence, only those metals that possess such 
an energy band may become superconducting at low tem- 
peratures. In fact, there is great evidence that this sup- 
position is true because I checked with 18 pure metals 
that superconductors and only superconductors possess 
a narrow, roughly half-filled a band in their (calculated) 
band structure.tH In particular, those metals (such as Li, 
Na, K, Rb, Cs, Ca Cu, Ag, and Au) which do not possess 
a narrow, partly filled a band do not become supercon- 
ducting. So far I have not found any exception to this 
rule. 

The modified BCS mechanism of Cooper pair forma- 
tion may also lead-to a better understanding of high-T c 
superconductivity.EJ In the one or two dimensional sub- 
lattices of the high-T c materials, phonons are not able to 
carry crystal-spin angular momentum. Hence, the nona- 
diabatic operator of C oulom b interaction in a a band, 
HqI, as given in Eq. (3.13) cannot depend on phonon 
operators. When a narrow, partly filled a band exists 
in such an anisotropic material, the electron spins are 
forced to couple to other crystal-spin- 1 excitations which 
must be sufficiently stable to transport the crystal-spin 
angular momenta. Likely, these excitations are coupled 
phonon-plasmon modes possessing a considerably higher 
energy than phonons. Therefore, the Debye temperature 
9 is considerably higher than in the isotropic materials of 
the standard superconductors and, hence, T c is higher. 
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APPENDIX A: SPIN DEPENDENT WANNIER 
FUNCTIONS 

The paramagnetic group of bcc niobium is given by 
M P = G + KG, (Al) 

where 

(A2) 



G = Ol 



is the bcc space group and K stands for the operator of 
time inversion. M p has the point group 



Go — Oh 



being the cubic point group Oh- 
The elements 

a = {a\R} 



(A4) 



(A5) 



of G consist of a point group operation a and a primitive 
translation R. The symmetry operators P(a) and the op- 
erator K of time inversion act on a function of position, 
f(r), according to 



and 



P(a)f(r) = /(a-V 



Kf^ = to, 



L i?) 



(A6) 



(A7) 



respectively. 

Consider the band structure of Nb depicted in Fig. [I], 
in particular the a band denoted by the dotted line. It 
is characterized by the four representations 



r' , H' N 2 , and P 4 



(A8) 



of O^. The Bloch functions of this single band cannot 
be unitarily transformed into usual Wfs which both are 
best localized and belong to a representation of Oh since 
^25; ^25> an d Pa are not one-dimensional Jffl 

However, such a unitary transformation becomes pos- 
sible in a natural way when we account for the existence 
of the electron spin. Let us define "spin dependent Bloch 
functions" by 



+ h 



(A9) 



where ^p^(r) denotes the Bloch function with wave vector 
k of the a band. The two-dimensional matrix [f sm (k)] 
is, for each k, unitary and 



u s (t) = S s t 



(A10) 



stands for Pauli's spin function with the spin quantum 
number s = ±| and the spin coordinate t = ±~. A sym- 
metry operator S(a) of the three-dimensional rotation 
group 0(3) acts on u s (t) according toll3 



S(a)u s (t) = u s (a 1 t) 

= ^d s > s (o!)iv(t) for a £ 0(3), 



(All) 



M p = G + KG , 



(A3) 



where the matrix [d s ' s (o!)] is the representative of a in 
the two-dimensional double- valued representation Di/ 2 
of 0(3). The quantum number m = ±5 of the crystal 
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spin distinguishes between the two functions belonging 
to the same wave vector k. If we have 



(A12) 



the two functions 4>j: m (r, t) with m = ±3 are usual Bloch 
functions with the spins lying in +z and — z direction, 
respectively. Otherwise, the functions <j>^ (r , t) still are 
usual Bloch functions with antiparallcl spins which, how- 
ever, no longer ly in ±z direction. 

The doubled-valued representations of the a band are 
given by 



D1/2 x V 25 



r+ + r+ 

D 1/2 x H' 25 = H+ + H$ 
D 1/2 x P 4 = P 7 + P s , 
D 1/2 xN 2 = N+. 



(A13) 



Hence, at the points T, H, P, and N the spin dependent 
Bloch functions of the a band can be transformed in such 
a way that at each symmetry point two functions form 
basis functions of the double- valued representations 

T+, H+, P 7 , and N+, 

respectively, of band 2 in Table ||. The spin directions of 
these basis functions, however, are different in different 
points of symmetry. 

In a bands, i.e., in the bands given in Table |, the ma- 



trix [f S m(k)] in Eq. (A9) may be chosen in such a way 
that 

(1) the Bloch functions <p^ m (r,t) vary smoothly through 

the whole k space and 

(2) the spin dependent Wfs 



"mi?- 1 1) = -L X: e"*%Jr, t) (A14) 

k 

are symmetry-adapted to M p according to 

(a)(r, t\T , m) 



BZ 



for a £ G and 

K(r,t\f,m) = ±{r,t\f,-m), 
where I have used the abbreviations 
(f,t\T,m) =w m (r- T,t) 

and 

T' = af+ R. 



(A15) 



(A16) 



(A17) 



(A18) 



The symmetry operators P{a) now act on rand t accord- 
ing to 

P(a)(r,t\f,m) = (a^ 1 r~a^ 1 R,a^ 1 t\f,m), (A19) 



with the meaning of a _4_being given in Eq. ( All[) . We 
define the plus in Eq. ( A16| ) to belong to m = +\- an d 
the minus to m — —\. The two-dimensional ma trice s 
[dm'm(ct)] are equal to the matrices [d s ' s (a)] in Eq. ( All ) 
and the matrices [d(a)d m i m (a)] = ±l-[d m ' m (a)] form the 
representation E, given in the first column of Table | (i.e., 
r^" for the band considered in this paper). The smooth- 
ness of the Bloch functions 4>^ m (r, t) guarantees that the 
spin dependent Wfs are optimally localizable. 

The spin dependent Bloch functions may be calculated 
from the Wfs by means of 



BvK 



ik-T 



w m (f- T,t), 



(A20) 



where the sum runs over the Born-von Karman volume 
(BvK). 

It is essential that the matrix [/ sm (fc)] cannot be chosen 
independent of k. This follows from the very fact that, 
on the one hand, the Wfs (r, t\T,m) would be usual spin 
independent Wfs if [f S m(k)] would be independent of k, 
and, on the other hand, we cannot assign to the consid- 
ered a band best localized, symmetry adapted, and spin 
independent Wfs since the representations r' 25 , H 25 , and 
P4 are not one-dimensional. 

An important consequence of the k dependence of the 
matrix [f S m(k)] is that Eq. ( Al5| ) does not hold when a 
symmetry operator acts on the spin coordinate t alone. 
That means, this equation does not hold when the oper- 
ators P(a ) are replaced by the operators S(a) given in 

Eq. (gnh, 



+ 4 



S(a)(r,t\f,m) £ ^ c m , m {a){r,t\f,m') (A21) 



for at least one a € Go,0 with 



^ |cw TO (a)| 2 = 1. 



(A22) 



Equation ( A21 ) holds since, otherwise, the functions 
(f, t\T, m) would be usual spin independent Wfs. 

The spin dependent Wfs form a complete basis in the 
a band. Thus, the functions S(a)(r, t\T, m) may be writ- 
ten as linear combinations 

S(a)(r,t\f,m) = ]T (f',m'\S(a)\f,m)(r,t\f',m'), 

f',rri 

(A23) 







with 



0\ = baa* a* ■ ct ct c,=s 

fli Q2 Q2 Qi 



(A31) 



E \(f',m'\S(a)\f,', 



(A24) 



belonging in Eq. ( A30 ) to SOS 1 . Oi does not belong 
to H' Cb since Ri ^ Q±. Hence we have 



and at least two nonvanishing coefficients 
(Ti, mi\S(a)\T, m) and (T2,m2\S(a)\T,m) belonging to 
different lattice points T\ and T2. Again, this statement 
holds for at least one a £ Gq. 



Equations (A15), (A16), and (A23) also may be writ 



ten for the fermion operators ct in Eq. (EO), 



[H' Cb ,S]^0. 
Finally, we derive the important equation 

fU-k) = ±f- s ,- m (k) 



(A32) 



(A33) 



which is needed in Sec. |III B| . The plus sign in Eq. ( |A33| ) 

holds for m = s and the minus for to — — s. 



P(a)ct P 1 (a) — d(a) E < K n , m (a)ct l for a e G, Equation (A26) leads with 



with \d(a)\ = E"L_i \dm'm(a)\ 2 = 1, 



Kct R- 1 = ±ct , 

Tm T,-m' 



(A25) 



(A26) 



and 



S(a)cJ, 5(a)- 1 = £ (T',m'|5(a)|f,m) ! 



ct 

T'' 



(A27) 



with at least two nonvanishing summands 
(Ti,mi\S(a)\T,m) and (T2,m2\S(a)\T,m) belonging 
to different lattice points T\ and T2. 

Now we show that, as a consequence of this equa- 
tion (A27), the operator H' Cb given in Eq. ( |3.3| ) does not 
commute with all the spin operators S(a). Reme mbe r 
tha t al l the matrix elements of H' cb satisfy equation (2.4) 
or (2.5) and consider a fermion operator combination be- 
longing to H' Cb , say 



O = J J 



C fi C f 2 C 'T2 C T 1 : 



(A28) 



where we have dropped the index to since it does not 
matter here. Further assume Eq. (A27) to be true for 
the point group operation a and the sum to consist of 
exactly two summands. With the abbreviation S = S(a) 
we may write 



Sets' 1 

T 



a-c^ + b-c^ 



(A29) 



where a ^ 0, b ^ 0, and Q ^ R. With Eq. flA29D we 
obtain 

SOS^ 1 = Set S^Sct S^Sc* S^Sa* S' 1 

T\ T 2 lz J-l 

x{a*c^+b*c &2 ){a*c^+b*c^). (A30) 
Consider, e.g., the operator product 



ct 

Tm 



1 



BZ 



v ? 

T 



ik-T t 

e cU 

km 



to the equation 



where 



Kct R- 1 = vim) • c f - 

km —k.—m 



v(±-) = ±l. 



(A34) 



(A35) 



Equation (3.23) yields the two equations 

E f- S ,-m(-k)cl n 



C 'k 
— k, — r 



k, — s 



(A36) 



and 



KC \m K "- E fsmikHs) ■ 

S_ 2 

becauseB 

JTct R- 1 = vis) ■ c f - . 

ks —k,—s 



(A37) 



(A38) 



Substituting Eqs. (|A36|) and ( |A37[ ) in Eq. ( |A35[ ) we ob- 
tain equation ( |A33j ) . 

Within the NHM, the localized functions (r,t\T,m) 
are replaced by nonadiabatic localized functions, 

(r, t\T , m) -> (r, £, g |f , to, i/) , (A39) 

which are orthonormal according to 

(T , to', n|T, to, n) 

= E] / / (T\m\n\r,t,(f)(r,t,q \T,m,n)drdq 



— fo'Tt'T&m' 



rpirp^m' m- 



(A40) 
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The qua ntum number v = n labels the functions satisfy- 
ingEq. (JU|). 

As a consequence of the third postulate of the NHM, 
also th e non adiab atic l ocalized functions satisfy the equa- 
tions ( A15 ) and ( A16 ). However, the symmetry opera- 
tors P(a) now act on f, t, and on the new coordinate q 
according to 



P(a)(f,t,q\T,m,n) 



"^a-^a-y |f,™,n). (A41) 



Consequently, the equations (A25) and (A26) are sat 



isfied for the nonadiabatic fermion operators m 
Eq. (^9|) as well as for the adiabatic fermion operators 

Tm 
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